Structure and the Lamb-shift-like quantum splitting 
of the pseudo-zero-mode Landau levels in bilayer graphene 



(N 

o 

(N 

in 



43 

I 



i 

c 

o 
o 



> 

o 

m 
m 

in 
o 

(N 



X 



K. Shizuya 
Yukawa Institute for Theoretical Physics 
Kyoto University, Kyoto 606-8502, Japan 

In a magnetic field bilayer graphene supports an octet of zero-energy Landau levels with an extra 
twofold degeneracy in Landau orbitals n = and n — 1. It is shown that this orbital degeneracy 
is lifted due to Coulombic quantum fluctuations of the valence band (the Dirac sea); this is a 
quantum effect analogous to the Lamb shift in the hydrogen atom. A detailed study is made of how 
these zero-energy levels evolve, with filling, into a variety of pseudo-zero-mode Landau levels in the 
presence of possible spin and valley breaking and Coulomb interactions, and a comparison is made 
with experimental results. 

PACS numbers: 73.22.Pr,73.43.-f,75.25.Dk 



I. INTRODUCTION 

Graphene, 1-5 an atomic layer of graphite, supports as 
charge carriers Dirac fermions that lead to unique and 
promising electronic properties. Of particular interest 
recently are bilayers 6 ' 7 (and multilayers) of graphene, 
where an added layer degree of freedom combines with 
spin and valley degrees to make the physics of graphene 
far richer. In particular, bilayer graphene has the prop- 
erty that its band gap is externally controllable. 7-11 

Quantum phenomena associated with Dirac fermions 
become prominent in a magnetic field. In a magnetic 
field graphene supports a set of characteristic zero-energy 
Landau levels, whose emergence and degeneracy have a 
topological origin in the nonzero index of the Dirac op- 
erators, or in the chiral anomaly in 1+1 dimensions. 12 
Monolayer graphene has four such zero-energy levels ow- 
ing to the spin and valley degeneracy. 

Bilayer graphene supports an octet of such zero-energy 
levels, due to an extra twofold degeneracy 7 in Landau or- 
bitals n = and n = 1. This orbital degeneracy brings 
about a new realm of quantum phenomena, 13-18 such as 
orbital mixing and orbital-pseudospin waves. The eight- 
fold degeneracy of the zero-energy levels is partially or 
fully lifted in the presence of Zeeman coupling, interlayer 
bias and Coulomb interactions, and these levels evolve 
into a variety of pseudo-zero-mode levels, or broken- 
symmetry states, as discussed theoretically in the con- 
text of quantum-Hall ferromagnetism 13 and others. 19 ' 20 
Experimentally the full splitting of the pseudo-zero-mode 
levels has been observed 21-23,25 in bilayer graphene. 

A key feature of graphene is that graphene is an intrin- 
sically many-body system of electrons equipped with the 
valence band that acts as the Dirac sea. Quantum fluctu- 
ations of the Dirac sea are sizable and lead to such quan- 
tum phenomena as velocity renormalization, 26 screening 
of charge, 27 and nontrivial Coulombic corrections to cy- 
clotron resonance. 28-32 It will be important to ask how 
such quantum fluctuations affect the pseudo-zero-mode 
sector in bilayer graphene, in order to interpret experi- 
mental results on broken-symmetry states properly, 



The purpose of this paper is to study the structure 
and spectra of the pseudo-zero-mode sector in bilayer 
graphene. It is shown, in particular, that the orbital 
degeneracy of the zero-energy levels is lifted by Coulom- 
bic vacuum fluctuations, leading to an appreciable shift 
and splitting of the n = and n — 1 levels; this is a 
quantum effect analogous to the Lamb shift 33,34 in the 
hydrogen atom. This vacuum effect is correlated with the 
Coulomb interaction within the pseudo-zero-mode octet 
to yield eventually a particle-hole symmetric spectrum 
for this octet. Unexpectedly, a detailed analysis of the 
interlayer Coulomb interaction reveals negative capaci- 
tance energies for bilayer graphene, which suppress pos- 
sible valley rotations in the pseudo-zero-mode sector. We 
discuss how this sector is split in spin, valley and orbital, 
with filling, in the presence of Zeeman coupling, inter- 
layer bias and Coulomb interactions, and compare with 
experimental results. 

In Sec. II we briefly review some basic features of 
the pseudo-zero-mode octet in bilayer graphene, and in 
Sec. Ill show that the orbital degeneracy of this octet is 
lifted due to quantum fluctuations. In Sec. IV we discuss 
in a simplified setting how this quantum effect cooperates 
with the Coulomb exchange interaction to determine the 
spectra of the pseudo-zero-mode levels. In Sec. V we 
take a close look into the valley breaking due to the in- 
terlayer Coulomb interaction. In Sec. VI we examine the 
hierarchy of broken-symmetry states under general and 
practical conditions with spin, valley and orbital break- 
ings, and in Sec. VII compare with experimental results. 
Section VIII is devoted to a summary and discussion. 



II. BILAYER GRAPHENE 

Bilayer graphene consists of two coupled honeycomb 
lattices of carbon atoms in Bernal A'B stacking. The 
electrons in it are described by four-component spinor 
fields on the four inequivalent sites (A, B) and (A',B') 
in the bottom and top layers, and their low-energy fea- 
tures are governed by the two inequivalent Fermi points 
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K and K' in the Brillouin zone. The intralayer cou- 
pling 70 = "fAB ~ 3 eV is related to the Fermi velocity 
v = (V^/2)a L j /h ~ 10 6 m/s (with a L = 0.246nm) in 
monolayer graphene. Inter layer hopping via the (A',B) 
dimer coupling 35 71 = -fA'B ~ 0.4 cV modifies the in- 
tralayer linear spectra to yield quasi-parabolic spectra 7 
in the low-energy branches |e| < 71. 

The effective Hamiltonian with such leading intra- and 
inter-layer couplings is written as 7 



where 



H 



bi 
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with p — p x + ip y , p' = Px — iPy Here = 
(ipA,ipB r ,i>A' iV'b)* stands for the electron field at the 
if valley, with A and £? referring to the associated sub- 
lattices; u stands for the interlayer bias, which opens a 
tunable gap 9 ' 10 between the K and K' valleys. We ig- 
nore the effect of trigonal warping oc 73 = jab' ~ 0.1 eV 
which, in a strong magnetic field, causes only a negligi- 
bly small level shift 27 ; also ignored are some nonlead- 
ing intra- and inter-layer couplings that lead to weak 
electron-hole asymmetry 36 in bilayer graphene. Hk is 
diagonal in the (suppressed) electron spin. 

The Hamiltonian Hk 1 at another valley is given by Hk 
with (v, u) — > (—v, —u), and acts on a spinor of the form 
\$> K — {ipB',i>A, i>B, ipA'Y- Actually, Hk 1 is unitarily 
equivalent to Hk with the sign of u reversed, 



H 



K' 



= s f n 



K 



(2.2) 



with S = diag(l, 1, — 1, — 1). In what follows we adopt 
Hk\-u for Hk 1 and simply pass to the K' valley by re- 
versing the sign of u in the if -valley expressions. Nonzero 
bias it^O thus acts as a valley-symmetry breaking. 

Let us place bilayer graphene in a strong uniform mag- 
netic field B z = —B < normal to the sample plane; we 
set, in Hk, P — > n = p + eA with A = A x + iA y = By, 
and denote the the magnetic length as I — \j\feB. It is 
easily seen that the eigenmodes of Hk have the structure 

*„= (|n)6W,|n-2) b%\ \n-l) b£\ \n-l)b^ (2.3) 

with n = 0, 1, 2, • • •, where only the orbital eigenmodes 
are shown using the standard harmonic-oscillator basis 
{|n)} (with the understanding that |n) = for n < 0). 

The coefficients b„ = (bP , b ( n ] , b£ } , b^f for n = 2, 3, . . . 
are given by the eigenvectors of the reduced Hamiltonian 
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(2.4) 



u c = V2v/e w 36.3 x w[10 6 m/s] t/B[T] meV, (2.5) 

with v measured in units of 10 6 m/s and B in tesla, is the 
characteristic cyclotron energy for monolayer graphene; 
M = 2 u/u c and 7' = ji/u> c . 

The energy eigenvalues obey the secular equation 

( 7 ') 2 (e' 2 - M 2 ) 

= {(e'- Mf - \n\] {(c' + Mf - \n\ + 1}, (2.6) 



where e' 



Let us denote the four solutions of 



the secular equation as eZ n < e_ n < < e„ < e+ for 
each integer n > 2 and \M\ < 1, so that the index ±n 
reflects the sign of e n . The eigenvectors b„ for \n\ > 2 
are written as 



{VH^l/(e' + M)}b^, 
- ( l/ 7 '){|n|-(e'-M) 2 }/ 
(e' - M)/VH 



(2.7) 



with fixed by normalization |b„| = 1. These expres- 
sions are equally valid for e„ — > ■ 

Of our particular concern are the n — and n = 1 
modes. For n = 0, "H re d has an obvious eigenvalue eo = 
\u = uj c M with eigenvector b = (1, 0, 0, 0)*. For n = 1 
Eq. (2.6) has three solutions (el-^ei,^ ) (excluding e' = 
— M). Actually le^J > 71 and we focus on t\ which 
is close to e . Let us set e\ — \u (1 — z), with z < 1 
determined from 



2 4(2- Z )(l-l¥)^ 



(7') 



(7') 2 + r 



(2.8) 



The associated eigenvector bi is written as 



bi = a (l,0, -(1 - M 2 z 2 ),-zM 



4 = 1/|1 + — {1- M z z z y + z 2 M< 



1-^.(2.9) 



Both z and z c , defined above, are functions of (M 2 , 7 ' 2 ), 
and are thus common to the K and K' valleys. They 
coincide for M = 0, z = z c = 2/(j' 2 + 1); numerically, 
z « 0.181 and c\ w 0.909 for M = at 5=10 T, with 
DRil.lx 10 6 m/s and 71 w 0.4 eV. These eo and e\ modes 
are called the pseudo-zero-modes since they evolve from 
the zero-energy modes of the M = case. For M = 
there arc eight such zero-energy Landau levels differing 
in spin, valley and orbital [n = (0, 1)] degrees of freedom; 
the presence of the zero-energy modes is dictated by the 
nonzero index 12,27 of the Dirac Hamiltonian Hk © Hk 1 - 
One can pass to the K' valley by setting M — > — M 
in the K valley expressions. The spectra (e ,£i) thereby 
change sign but for later convenience we continue to use 
n = (0,1) to specify the pseudo-zero- mode levels at the 
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K' valley; we thus write e n \K> = —£ n \K for n £ (0,1). 
The eigensystems at the two valleys are related as 



H hl projected to the octet of pseudo-zero-mode Landau 
levels is thereby written as 



e„k' = -e-nk, (2.10) 

for \n\ > 2 [and for n £ (0, 1) as well if one sets ±0^0 
and il4l]. 

The Landau-level structure is made explicit by passing 
to the \n, yo) basis (with yo = £ 2 p x ) via the expansion 
(**(x),**'(x)) - E„,, (x|n,y ){C ;a (2/o)}, where n 
refers to the Landau level index, a £ (t, I) to the electron 
spin and a € (if, if') to the valley. The charge density 
P-p = I d 2 xe J P x p with p = (tf*)ttfX + (tf*")t$*" i s 
thereby written as 27 



P-p - 7 P 2^ Z^ 5 p 



oa 
P - 



R 



kn;ab 
a/3;p 



k.n— — oo a, a 



| dy ^ ;at (yo) 



)e^ r ^ ;h (2/o), (2.11) 



where 7 P = e~^ 2p2 / 4 ; r = {il 2 d/dyo,yo) stands for the 
center coordinate with uncertainty [r x , r y ] — it 2 . In par- 
ticular, the charge operators ^ a i?Q™' ; p° obey the 
algebras 38 associated with intralevel center-motion and 
intcrlevel mixing of electrons. 

The coefficient matrix g^ l ' a = g^ l \ a at valley a £ 
(K, K') is constructed from the knowledge of the eigen- 
vectors b n | OJ 



9 k p n 



where 



6«6W/W.N+6«6W/ 1 !*l- 2 .N- a 
HbFbW+bPWfW-W- 1 , (2.12) 

^^(^)""^(^ 2 ) (2-13) 

for k > n > 0, and f p lk = (/*™)t; p = p x + ip y ; it is 
understood that = for fc < or n < 0. As seen 
from Eq. (2.10), they are related at the two valleys so 
that 



mn;K f — m,— n\K 



„mn:a | 

5 P |m 



-Mi 



(2.14) 



where it is understood that one sets ±1 — > 1 and ±0 — > 0. 
Within the n e (0, 1) sector, they are common to both 
valleys, 



1, ffp 1 



1 



„2 l/>2„2 

-"1 2 



^ 2 p 2 , 



,10 



£p/V2, (2.15) 



so that the valley index a may be suppressed. 

From now on we frequently suppress summations over 
levels n, spins a and valleys a, with the convention that 
the sum is taken over repeated indices. The Hamiltonian 



H u — e™ 8R 



nn 



Pi (T 3 ) l3a R i 



a/3;0 



(2.16) 



with n £ (0,1), where 6° = \ 

fift-mri — j^rnn;KK v mn;K' K' 



R. 



t and t\ — \ u (1 — z); 

H o = « a( 5. ~~ J - L a p-o ■ Here the Zeeman term 
pz = g* /j,qB is introduced via the spin matrix T 3 = 
As interlayer bias u is turned on, these n = (0,1) levels go 
up or down oppositely at the two valleys. Nonzero u thus 
critically breaks the valley symmetry of the n = (0, 1) 
sector. 



III. VACUUM FLUCTUATIONS 



The Coulomb interaction is written as 
V = ^J2 v p : P-pPp 



(3.1) 



where t> p = 27ra/(eb|p|) is the Coulomb potential with 
a = e 2 /(47re ) ~ 1/137 and the substrate dielectric con- 
stant £b5 J2 P — I d 2 p/ (27r) 2 . In this section we show that 
the orbital degeneracy of the pseudo-zero-mode octet is 
lifted by Coulombic quantum fluctuations of the valence 
band (the Dirac sea). 

Let us define the Dirac sea |DS) as the valence band 
with levels with n < — 2 arc all filled. We take the expec- 
tation value (DS|y|DS) to construct the effective Hamil- 
tonian that governs the electron states over |DS). The 
best way to achieve this is to construct the Hartree-Fock 
(HF) Hamiltonian V KF out of V. 

In this paper we generally focus on many-body ground 
states \G) with a homogeneous density, realized at in- 
teger filling factor v £ [—4,4], and set the expectation 
values (G\R™pj£ b \G) = <5 k , Po v^ ab with p = l/(2n£ 2 ) 
and <5k, o = (2tt) 2 <5 2 (k). Accordingly, the filling factor 
v nn\aa _ ^ w hcn the Landau level specified by (n, a, a) 
is filled up. 

Let us divide V HF into two parts, V HF = V D + V X - 
The direct interaction reads 



Vb = PQVp^o v\ 



m'm';bb 



nn;aa r> 



(3.2) 



while the exchange interaction is written as 

t/ \ ^ 2 mn'-,b m'n:a , mn;ba r/m'n';ab. /,, o\ 



note that — 5 mn . In Vb and Vx we sum over 

filled levels (to, n) and retain the pseudo-zero-mode sec- 
tor m',n' £ (0, 1). As usual, the direct term Vq with a 
divergent factor w p ^o is removed if one takes into account 
neutralizing uniform positive background charges. 

Let us first extract, out of Vx, the contribution from 
the Dirac sea, i.e., all filled levels with n < —2, 



yUS 



I>p7p 

p 



E 

n<-2 



\ n m'n;a\2 rym'm':aa /q a\ 

\9p I H aa;0 ) ( 6A ) 
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where the sum over m' € (0,1), a e (K,K') and a <E 
(t,-i) is understood. Actually, the sum over infinitely 
many filled levels with — Nl < n < — 2 gives rise to an 
ultraviolet divergence with Nl — > oo. 

Fortunately one can isolate the divergence and even 
evaluate Vj? exactly for M — > 0, as shown below. It 
is clear from Eq. (2.14) that g™ n '> K = g™ n ' K ' = g~ m ~ n 
for M — 0. For M — > one can thus replace the sum 

E„<- 2 -> 3 (E„<-2 + E„> 2 ) in E q- ( 3 - 4 ) and note thc 
completeness relation 39 

oo 

E ISp "I 2 = e £2p2/2 (3.5) 

n— — oo 

to obtain 

E l3p"| 2 M =° ^(^ p2/2 -|5p°| 2 -|5 p 1 | 2 ), 

n<-2 

E ICI 2 M =° ^(^ 2p2/2 -|5p°| 2 -|5 P 1 | 2 ). (3-6) 

n<-2 

We have independently confirmed this result by direct 
numerical calculations. 

The e e p / 2 term in Eq. (3.6) leads to a divergence 
upon integration over p; it is, however, common to all 
levels m' and is safely omitted. We thus take — \ (|g p °| 2 + 
l^p 1 ! 2 ) as the regularized expression for V^„<_ 2 l.9p"| 2 ; 
and -\ (l^l 2 + | 3p °| 2 ) for £„<_ 2 l^pT- One"can make 
use of these M = expressions to evaluate the M^O 
corrections unambiguously. 

Substituting these regularized expressions into 
Eq. (3.4) and integrating over p yield 

x/DS o p00;aa 1 pll;aa /n 7 x 

- f vV;0+ £ vVo. W-l) 

with 

e ° = E^pKI5p i 2 + i< 1 i 2 ) = H 1 + H)^- 
p 

4 = E^pKl5 P 1 | 2 + l5 P °l 2 ) = l( 1 + I c2 -^)^ 
p 

C = c 2 (4-3c 2 ), (3.8) 
where 

E«p7p = yfr/2V c = V c , 
p 

K = a/(e b £) w (56.1/et) /B[T] meV. (3.9) 

Numerically, for 71 = 0.4 cV, u = 1.1 x 10 6 m/s and 
M = 0, one has c\ = 0.909 and C = 1.156 at B = 10T, 
which give e° w 0.73 V c and 4 w 0.58 V c . 

Vacuum fluctuations from the Dirac sea thus shift and 
split thc n = and n = 1 levels; this is somewhat con- 
trary to the common wisdom that filled levels can be 
regarded as dormant, with no active role to play. The 
energy splitting 



reflects the difference in their spatial distributions, as is 
clear from Eq. (3.8). Thc n = levels become higher in 
energy than the n = 1 levels, when they are unoccupied. 
Actually their relative positions vary with the filling of 
the n = (0, 1) sector. To get a rough idea about how 
they vary, let us suppose filling the n = 1 levels first (for 
u = (iz = 0). Equation (3.4) then tells us to include 
extra contributions — |g 01 | 2 and — I 1 1 1 2 for e° and e*, 
respectively; this yields e° = \ (1 — \c\) V c and el — 
~\ (1 — l c i — 3^) ^ ms tead, the n = levels were 
first filled, one would find e° = — \ (1 — \c{) V c and = 
2 (1 — \c\ — jC) V c . The ground-state energy is thus lower 
in the latter case. This shows that the level structure 
indeed varies with filling of the sector and, in addition, 
indicates the unusual character of the orbital-polarized 
ground states, which we examine in detail later. 

Let us next suppose that a pair of n = and n = 1 
levels of the same spin or valley were filled (in the pres- 
ence of some spin or valley breaking). One would then 
find that e° and e* reverse sign so that the filled n = 
level is now lower than the filled n = 1 level; this leads 
to a gap 2el = (1 + ^c 2 — jC) V c between the filled and 
empty n = 1 levels. This (spin or valley) gap is larger 
than the possible "orbital" gap (1 — \c\ — \G) V c between 
the n = 1 and n — levels, discussed above. This sug- 
gests that relatively large energy gaps expected for the 
v = ±2 and v = ground states are either spin or valley 
gaps rather than orbital gaps, in accordance with Hund's 
rule and Ref. 13. 



IV. HF CALCULATION 

In this section we examine the level structure of the 
pscudo-zero-mode sector more quantitatively. To this 
end, let us extract from Vx in Eq. (3.3) the exchange 
interaction acting within the n = (0, 1) sector, 

vF = -E^p°' 
p 

O r I ^00(2 00;6a , i 10|2 „ll;6al p00;a6 

U - llffp I V f) a +\9p\ Vp a )^ a /3;0 

J-/U 10 l 2 , 00;ba 1 11 12 , ll;6o-i p ll;o6 
+ ilSp I v fi a + \9p I Vf) a )^a/3;0 

_L^°° r, 11 /,.10;6a p01;ob , 01;6o plOjabi r A -\\ 

We consider H u +Vj? S +V£ z , as the effective Hamiltonian, 
to explore the pseudo-zero-mode sector. The exchange 
interaction Vj? S + V£ z is invariant under rotations in val- 
ley and spin, but not in orbital (0, 1) space. In contrast, 
H u , consisting of inter layer bias u and the Zeeman energy 
/iz, can lift all of the valley, spin and orbital degeneracies. 

The vacuum orbital splitting Ae v = \CV C , u and /iz 
would in general compete to determine how the pseudo- 
zero-mode octet is split in orbital, valley and spin when it 
is gradually filled with electrons. In practice, a relatively 
weak bias u can easily overtake thc Zeeman energy 



Ae v = e° v -el = \CV C 



(3.10) 



m = 9*VbB w 0.12B[T]meV = 1.4E[T]K; (4.2) 
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(i) fiz > (l-z)u 



(ii) (l-z)u > Hz 



0+1 



t U +T 



0_t 







l_t +T 



l_l 1+1 



1+T 



-T 



FIG. 1: Empty levels in the pseudo-zero-mode sector at v = 
—4; for illustration, (1 — z)u/hz = 0.2 for (i) and 1.5 for 
(ii), with h^/Vc = 0.05. The dotted lines refer to the energy 
scale £° = \ (1 + | cf) V c . The valley index ± here stands for 
(+,-) = (K,K') with u > 0. 



we allow u to be arbitrary but keep \u\ < V c . See Fig. 1, 
which depicts the empty n = (0, 1) sector (at v = —4) 
governed by H u + V^ s , with level spectra 



n± _ n i r, 
t i — fc v =n fc„ 



n± n _i_ F r. 

t -|- — fc V ^ C l 



(4.3) 



where e" a stands for the energy (per particle) of the n aa 
level with orbital index n £ (0, 1), valley a £ (+, — ) and 
spin a £ (t,+)- For valley indices we use (+, — ) to spec- 
ify that the (-) state is lower in energy than the (+) state 
for each n and spin a; accordingly, (+, — ) = (K,K') for 
u > 0. There are two possible level patterns, depending 
on (i) fiz > (1 — z)\u\ or (ii) (1 — z)\u\ > Hz- For defi- 
niteness, in this section we assume Hz 3> M or even set 
u — > 0, and examine the hierarchy of broken-symmetry 
states of case (i); we discuss the general case with both 
u and Hz later. In this setting, H u practically has little 
effect on the ground-state energies (since Hz/V c <C 1). 
Accordingly, and for further simplification, we keep in 
mind that H u triggers spin breaking, suppress it and fo- 
cus only on V = V® + in what follows. 

To diagonalize the exchange interaction V we proceed 
as follows: Note again that V is invariant under rota- 
tions in spin and valley. Spins and valleys, if driven ex- 
ternally, can therefore easily rotate. Obviously, under 
the magnetic field B z = —B, spins remain polarized in 
(t,+). Valleys may, however, rotate, (K,K') — > (+,—). 
We suppose that the spin (t, +) and the rotated valley 
(+, — ) are good quantum numbers, and diagonalize each 
(valley, spin) sector via rotations in orbital n — (0, 1) 
space. We thus rotate ip™' a = (V>a a > V'a") nrst in valley 
space and then in orbital space, 

*PT a (yo) = [U(0, 4>)]\ [U(d, 0)]™ ^\y ), (4.4) 
using the SU(2) matrix 



cos(0/2) -e-^sin(0/2) A 
e^sin(0/2) cos(0/2) J 



(4.5) 



By construction the orbital rotations (0, 0) refer to each 
(valley, spin) sector, i.e., (0,0) — > (9 aa , 4>aa)'i f° r con- 
ciseness, we suppress this reference. For the transformed 



fields <£>™ ;a we use n = (0, 1) and a = (+, — ). For exam- 
ple, $° ;+ stands for ip^ K for 9 = = while it stands 
for ip^j K for = 7r and = 0. Actually, the valley rota- 
tions (K, K') — s- (+, — ) are left undetermined for V alone. 
They are to be fixed later (in Sec. V) when we consider 
H u and a valley breaking due to interlayer Coulomb in- 
teractions. Note next that and can be absorbed into 
the relative phases between $^ ;a and , we therefore 
set = ^ = from now on. 

T 0;ba /^<i^0;bt- 



Let N' 



/3ct 



(G\^ b %^\G) and N 1 * 



■x 



(G|$£ 6t $^ ;a |G) denote the filling factors in the form of 
matrices in valley and spin. By our assumption they are 
diagonal in both valley a £ (+, — ) and spin a £ (t, i). 
Via rotations (4.4), Vj? S reads 

^ X DS = e°(0) TC^ + 4(9) + l^m' X Z, (4-6) 



with 



4(0) = ^[±(l + ±c 2 )-iC(l-cos0)], 
4(0) = y c [I(l + I c 2)_ 1 Lc(l + cos0)]. (4.7) 



b denote the charge operators for ^' b , i.e., 



iC'-t with ^ b -> and x$ = n\ ,. r 

[. . .]' = {d/d9)[- ■ •]. Similarly, V£" reads 



01;afc , -t-j 10;a6 



+ [E* (0) N°* a + E x v (0) N$ a ] K 1 ^ 



+ 



1 ^(0) N°* a - El, (6) N^ a ] ' X%, (4.8) 



where 



£;°(0) = -f c [l-^C(l-cos0) 2 ], 
E° v (0) = ^(0) = -K[icf- 1 LC , sin 2 0], 
E v {6) = -V4[l-^C(l + cos0) 2 ]. (4.9) 

See Appendix A for details. It is clear now that and 
V£ z are divided into four sectors, specified by (a, a), of 
the 2x2 matrix Hamiltonians for ($° ;a , $^ a )- 

Let us start filling the empty n = (0, 1) sector at 
v = —4, which consists of four n = and and four 
n = 1 levels of energy e°(0) = \ (1 + \ c\) V c and 4(0) = 
e°(0)-±CK; numerically, (eJ(Q),e°(0)) = (0.58, 0.73) V c 
at B = 10T. See Fig. 1. There the orbital breaking 
Ae v singles out the $j~ level (or n aa = as the 

lowest-lying one, which will thus be filled first. Sup- 
pose that this level is filled with fraction n\ < 1 and 



substitute iVi! 



= ni and 



into V. The 



(—,|) sector then consists of {el(9) + niE v (8)}TZ,,'. Q , 

{e°(0) + 7ii£^(0)}7?.°°: o ~ and the (0, l)-mixed compo- 
nent 



-[£(9) + % ni fy(9)]'x 



n 



(4.10) 



Accordingly, on choosing so that this term dis- 



appears, the (— ,i) sector is diagonalized. Solving for 



6 




tor with ($ 1; ~, <f> 0:_ )|e=7r = , of energy 

(— 4(0), — ej(0)), and is polarized in valley and spin, with 
a gap 24(0) [~ 1.17 K for S = 10T]. 

Among the empty levels at v = — 2 the lowest one is 
the 1+4, level, rather than for fiz > U. To reach the 
v = — 1 state, one therefore has to fill this l+j, level. The 
analysis is essentially the same as for the v — — 3 case, 
and one finds the filled l m +j, and empty m +4, levels of 
energy =p~ (1 — ~ c 2 ) V c . The v = — 1 state is polarized in 
valley and orbital, with an orbital gap (1 — | c\) V c . 

One can further go to the v = state by filling the 
m +i level, which thereby turns into the filled 0+4|e =7r = 
1+4, level. The v = state is polarized in spin, with a 
gap 2 4(0). 

Repeating essentially the same analysis for the spin-up 
sector takes one to < v < 4. The resulting spectrum is 
electron- hole symmetric, as shown in Fig. 2 (b). 



FIG. 2: (a) Variations of the spectra of the <f>{ and $° levels 
in the range v £ [—4,-2]. Note the symmetry e $ o- (9) = 
— e $ i- (tv — 9). (b) Spectra of the pseudo-zero-mode Landau 
levels at each integer filling factor v 6 [—4,4], with pz > 
u > and pz +0. Blue blobs refer to spin-down levels 
(1-4., 0-j,, 1+4., 0+4,)^ from left to right and red squares refer 
to (1-f, 0-f, 1+f, 0+t))|e; empty levels have positive energy 
and occupied levels have negative energy. 



V. INTERLAYER COULOMB INTERACTION 

The Coulomb potential between the two layers slightly 
differs from the intralayer potential u p , and leads to a 
weak valley-SU(2) breaking. In this section we examine 
the effect of this breaking. 

The relevant Coulomb interaction is written as 



then reveals that 6 = for < n\ < 1/2 while 6 rises 
to 7r/2 as 7ii is increased from 1/2 to 1; see Fig. 2 (a). 
We thus find that the filled $}~ and empty $?~ levels at 



-I 



v = — 3 are equal mixtures oc tp, 
and n = 1 orbitals, to be denoted as 1 



of the n = 



-4- 



1_ 



and m _4 = 1 =7r /2 ? respectively, in obvious notation. 
The v = —3 ground state thus consists of the filled l m _4 



and empty m „4, levels of energy ±i (1 



2 °i 



V c , with 



an orbital gap (1 



2 °1 



) V c [« 0.55 V c for 5 = f 0T] 



p 



(5.1) 



where = i v p (1 ± e rf ' p ') with the interlayer separa- 



J P 2 

tion d w 3.34A. Here p~ = p uppor - p lowcr stands for the 
charge-density difference between the upper and lower 
layers. One can write it as p~ = p~' K ~ p~' K with p~' ,a = 
^ T ya for a e ^ ^ and r = diag(l, -1, -1, 1). Let 



us denote their Fourier transforms p_p = J cf x ( 
as 

IP J-P=7p 2^ K aa;p > etc -: 



ipx 



/0" 



(5.2) 



As seen from Fig. 2 (a), the 0_+ 1 e level has come down 
along with I-4I0, owing to exchange interactions. One 
can therefore reach the v = —2 state by filling the 0-i\g 

level, i.e., by setting A^|"" = 1 and AT°j = n with where h* n < K = h* n \ K are given by g* n \ K in Eq. (2.12) 
< no < 1. Diagonalizing the (— ,+) sector then shows 
that 9 varies from 7r/2 to 7r as uq is increased from 
to 1/2, and stays at ir until no reaches 1. Note that, 
as one passes from v = — 4—!-— 3 — ^ — 2, the lowest 
energy level <fr^'~ has undergone the change 1_ — > l m _(= 
l_| e=7r / 2 ) — * 0_(= l_|e =7r ). Here we notice the rule 



with the signs of the and &!? &r?^ terms reversed. 



Within the (0, 1) sector h% 
^ P 9p 



fen 



hold, except for 



3 — n „2 



that each empty level characterized by $™ ;± |e = o = Va ;± 
turns into ^"^I^Tr when it is filled up. Note, in this 
connection, the relation 

4(0)+£^(0)+^(0) = -4(0) (4.11) 



2cf 



(1 - M 2 z 2 ) 2 



7 



12 



(5.3) 

= z c -2c\z 2 M 2 - (5.4) 



and an analogous one for £v(^)j which imply 
that the filled and <&^~ levels have energies 

(-4(^,-4(70) = (-4(0), -4(0)). The v = 2 
ground state thus consists of the filled (—,4) sec- 



it is thus unnecessary to specify their valley indices. 

Let us now substitute w+ = v p — v~ into V, extract 
the 0(v~) correction AV to V, and construct the HF 
Hamiltonian from it. We first consider the Dirac-sea con- 
tribution AV DS 
reads 



A1/ D DS - 



AlAp s . The direct interaction 



AV^ S = W 



n<-2 



5h^ n 



m=0,l 



(5.5) 



7 



where <5/iq™ 



unn I unn I 



which is odd in M cx 



1//9 Qnrl AB mm — r>mm;KK r>mm;K' K' . 
U/Z, dllU OHpp. Q — Mpp. rt /90;O 



w ee y c » 0.0104 v^PTv; 



(5.6) 



comes from /9ow p= o = 5 W. From AV^ we have omitted 

~,mm\bb 



a term cx WR 



/3/3;0 



that shifts all levels uniformly. 



This AV^ b leads to valley polarization cx hg m SR^ 
of 0(MW) in the n = (0, 1) sector. A direct numerical 
calculation shows that 



E 

n<-2 



5h™ 



-16.0 M 



(5.7) 



for M < 1 at B = 10T. AVg s thus acts like H u 
and, when combined with it, effectively reduces u — > 
(1 — 16W/w c )u. With this in mind, we henceforth un- 
derstand that H u already involves the effect of AVj!p s . 

The exchange interaction , on the other hand, 

takes the form of V® s in Eq. (3.4) with replacement v p — > 



~,m'n:a 1 2 



This AV? S is 



v- and "i"^ -+ \h% n -' a \ 2 -\g p > | . ^ ^ x 
calculable numerically. Fortunately, for M = 0, one can 
evaluate it exactly, as done for V^ s , with the result 



DS M- 



X 



>o 1 \ ^ 



- \9 P 



2 )i& (5-8) 



which, upon integration over p, turns out to vanish to 
0(W). We have also verified this result numerically. 

The main 0(W) corrections come from the interaction 
AV pz acting within the (0, 1) sector. See Appendix B for 
the details. Here we quote only the result : 

AV pz = IWi-vU™** + I{6) + t{8)zl 1 } 

+\W [K(9) + t'{9) zl 2 ] - t{6) T w + ■ ■ ■ , 
IV = {e: ap (0)SN°p a + ^(6)6^^611^ 
+\ {e^{9) SNp a - 4 p ((9) 5X a p 
+W{AI(9) + \AK(6)}, (5.9) 

with n e (0, 1), t{6) = ±(sin6») 2 and 

1(9) = [r 5N° a + ri SNi a ] [r Q 5n%. + nM^o] , 
K{9) = [{rl)'5Nl a -(rl)'5Nl a ]5Xpp 

X 2 = (6 01 N++) X# + (5^N--) X-p~, (5.10) 

where ro ee r$(9) = 1 — zsg and n ee r\(6) = l — zcg with 
se = sin(6>/2) and c e = cos(6»/2); r' = (d/d9)r (9), etc. 

Here S denotes the difference between the (++) and ( ) 

components so that 5N% p = AT°j++ - JVj,", 5ftJ$| ee 

- ^b" = - ^7' ctc - Similarly, 

<5 01 denotes the (0, 1) difference: S 01 N++ = iV°° ;++ - 

^Kto = <°^o + -^o + ' etc - and 
AK(9) stand for the crossed-spin variants of 7(0) and 



K (9) such that Atf(0) = [■ • -] n 5X u + [■ • -J^Tt. etc - 
From AF pz we have omitted terms (• • •) cx t'(9) that 
eventually vanish when t(9) is chosen to be stationary. 

The valley-dependent term —t(9) Y\y in AV pz involves 
capacitance energies, defined as 



e° ap (0) = WD(l-co S 9) 2 , 

& P (*) - e° cap (9)^W[l-cl-z + D( S m9n 



= W^(H-cos0) 2 , 



(5.11) 



where D = ± cf (1 - c 2 ) - I* + ±z 2 -± z + |z 2 < 0. 
Here we encounter negative capacitance energies 



e" (0)<O, 6" (0)<O. 



(5.12) 



This has an important consequence: It is easily seen that 
(G|IV-|G) < for the possible ground states |G) dis- 
cussed in the previous section. The capacitance energy of 
the n = (0, 1) sector therefore is minimized for t(9) = 0, 
or 9 = 0, 7T, which implies that there arise essentially no 
valley rotations in the pseudo-zero-mode sector in bilayer 
graphene. This is in sharp contrast to conventional quan- 
tum Hall (QH) systems, for which positive capacitance 
energies favor valley-symmetric and antisymmetric com- 
binations (9 = it/ 2) for eigenstates. Normally both direct 
and exchange interactions contribute to the capacitance 
energy and their contributions tend to cancel for zero 
layer separation d ~ > or to 0(W); this leaves positive 
0(d 2 /£ 2 ) capacitance energies for conventional systems. 

These negative capacitance energies for the n = (0, 1) 
sector are traced back to the fact that the n = 1 modes, 
in part icula r, are distributed on both layers with ratio 
1 to \J 1/2 in amplitude, or with ratio 1 — z/2 to z/2 
in charge, as seen from Eqs. (2.9) and (5.3). The di- 
rect Coulomb interaction thereby gets somewhat weaker 
(ro,ri < 1) and is overtaken by the exchange contribu- 
tion, leaving (e™ap'^cap) negative to 0(W). This fact, 
though clear in terms of 4-component electron spinors 
^ K and ^ K , could easily be overlooked when one uses 
the approximate 2-component spinor formalism 7 for the 
description of bilayer graphene. 

We study the effect of this 0(W) valley breaking as 
a perturbation to the eigenstates found in Sec. IV. One 
may set t(9) = and read from Eq. (5.9) the 0(W) 
corrections to the spectrum: 



i± 

Wi 



p0± 



= \W{ 
= \W{ 



i/iron 5N° aa ±rl 5N^ a 
± r r 1 SN^ a 



v±r 2 8Nl 



zt{9)8^N^}, 
+ zt(e)6^N±±}; 



'44 
(5.13) 



analogously for the spin- up sector. Here E 



stands for the correction to the 1£ 



U;0 



sector. For the 



v = — 3 state with only one filled level, e.g., one finds 



1 _ 0=7T/2 

Wi = 

, _ 9=k/2 
Wi " 



(2 



1 =2 



z 2 )W, 



8 



E 



1+ 



9=0 



= (l + z-±5 2 )IF, 



VI. GENERAL CASE 



(5.14) 



We have so far supposed that fiz ^> u ~ 0. Let us now 
relax it and consider the full effective Hamiltonian 



H cS = V + H u + AV, 



with V = V£ s 



Vl z and AV 



AV» S 



AV» S 



(6.1) 



AV pz 



to explore the pseudo-zero-mode sector with both fiz and 
bias u. We keep |u| <C V c so that we can still use the u = 
expressions for V + AV, with u retained only in H u as a 
small perturbation. [We have in mind u ~ a few meV at 
B ~ 10 T. ] Actually this is a reasonable approximation 
if one notes the following: Possible 0(u) valley breakings 
that come from the Dirac-sea contribution (V® s + ■ ■ ■) 
can effectively be taken care of by rescaling u in H u , 
as we have seen for AVg s in Eq. (5.7). In addition, 
ypz _|_ AV pz , acting within the (0, 1) sector, contains no 
O(M) corrections; note in this connection that z c and z 
differ from z only to 0{z 2 M 2 ). 

The energy scale ~ V c of the n = (0, 1) sector is 
set by V, and relatively weak (spin, valley and orbital) 
breakings in H c g trigger symmetry breaking of this sec- 
tor, leading to exchange-enhanced gaps at each integer 
filling factor v e [—4,4]. We have seen that negative 
capacitance energies suppress valley rotations. Accord- 
ingly, one can assign (+,— ) = (K,K') or 9 — 0, as- 
suming, without loss of generality, u > 0. Via rota- 
tions (4.4), H u somewhat changes its form; the main 
change is e™ 5Bffi. -> e n u {9) 51V$. + ■■; with 



4(*) 



| UCOS0{1 
I UCOS0{1 



-cos 61)}, 
iz(l + cos6>)}. 



Let us look at Fig. 1 again, which depicts the empty 
n = (0, 1) sector for the two cases, (i) (1 — z)u < [iz and 
(ii) (1 — z)u > /iz; see also Fig. 3. There the lowest-lying 
level is in both cases and is filled first. As the 1_| 
level is being filled, it comes down in energy, followed 
by the 0_j^ level (of the same valley and spin) paired via 
the exchange interaction. The v = —3 state will therefore 
consist of the filled l m -j, and empty m _^ levels of energy 

^ = -i(l-i C ?)V c +[6i(0)+<-J e=V2 , 

= +W-H)V c + [el{e) + E^] e ^ /2 , (6.3) 
which lead to an orbital gap, 

efi-3 = (1 - -A) Vc - \zW. (6.4) 

As a result, the v — —2 state is also unique and consists 
of the filled (l_j., 0_^)|e =7T levels, polarized in both valley 
and spin. (If, instead, an empty n — 1 level were to 
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FIG. 3: Spectra of the pseudo-zero-mode sector for each in- 
teger filling factor v £ [—4,4], with fiz/V c ~ 0.027 taken as 
a typical value at B = 10 T; this yields u CI /V c ~ 0.093. (a) 
u = 0.1 ti cr . (b) u = 0.5 iter, (c) u = 1.2ti cr . Blobs and 
squares denote levels ordered as in Fig. 2 (b). 



come down in passing to v = 2, it would lead to an 
orbital-polarized v = —2 state, which is energetically not 
favored, as noted at the end of Sec. III. We therefore 
exclude this possibility.) 

Among the empty levels at v = — 2 the lowest- lying 
level is either 1_^ or 1 + ^, depending on case (i) and (ii). 
Comparing their energy levels reveals that is lower 
for u < u CI , and vice versa, with the critical value it cr 
given by 



(1 - z) u CI = /i z + X(z) W, 



(6.5) 



where X(z) 



35 



~2 



Cases (i) and (ii) are thus 



distinguished by this more precise criterion u < u CI and 
u > iter- Correspondingly, the v = — 2 level gap differs in 
structure for the two cases, 



e i/=-2l(i) 



2el(0) + (l-z)u-X(z)W, 
fiUlfii) = 2 4(0) +^z. (6.6) 



The v = — 2 gap changes from a valley gap to a spin gap 
as u is increased across u rr - 



9 



Numerically, at B = 10T, X(z)W/V c w 0.049 and 
Mz/Vc ~ 0.027 for typical values ^ z ~ 1.2 meV and 
£{, <~ 5; this leads to u cr ~ 0.093 V c ~ 4 meV; or more 
generally, u cr /V c - (0.065,0.13,0.17) for B = (5,20,30) 
T. Note that the Coulombic inter layer valley breaking 
X(z)W is practically comparable to hz- 

Filling cither or l_-f leads to the v — — 1 state. 
It consists of the filled l m +^ level (and the paired m+ jj 
[over the filled v = —2 basis (l_j,, 0_j_)|e =7r ] for u < u cr 
while it consists of the filled l m _f level (and the paired 
m -t) for u > u cr . The resulting level gap ef, a ?_ 1 , never- 
theless, turns out to be the same as 

The v — state is reached by filling the paired n = 1 
level, and differs in composition, depending on u; see 
Figs. 3 (b) and 3 (c). It is spin-polarized for u < u CI 
while it is valley-polarized for u > w cr , with a level gap 

ef = p | (i) = 2 4(0) + /i Z -(l-*)u, 

4 a = P ol(ii) = 2el(0) + (l-z)u-fi z -2X(z)W.(6.7) 

It will be clear now how to reach the v = 1 <~ 4 states. 
One eventually finds that the pseudo-zero-mode octet has 
perfectly particle-hole symmetric spectra, as shown in 
Fig. 3 for each integer filling factor v e [—4,4]. It is 
evident from the figures how each level behaves as the 
n = (0, 1) sector is gradually filled, and that the spectra 
change in pattern for u < u cr and for u > u cr . 

It turns out that the level gaps at odd integer filling are 
all purely orbital gaps of the same magnitude, e52±3 = 
e^?±i, which are relatively small and insensitive to both 
Hz and u. It is interesting to examine experimentally how 
these possible gaps respond to a tilted magnetic field. An 
additional parallel magnetic field B\\ effectively works to 
increase Hz and u CI . Equation (6.4) therefore suggests 
that the v — ±1, ±3 gaps are insensitive to Bn. 

So far we have supposed filling the n = (0, 1) sector 
gradually for a fixed u. When u is applied for a fixed 
density or v, cases (i) and (ii) are simply connected for 
the v — ±3, ±2 states but not for the v = 0, ±1 states. 
We discuss this point in the next section. 



VII. COMPARISON WITH EXPERIMENTS 

Experimentally full splitting of the pseudo-zero-mode 
Landau levels has been observed 21-25 at high magnetic 
fields. The quantum numbers (especially, valley and or- 
bital ones) of the resulting broken-symmetry states re- 
main unknown yet, but the way they emerge with mag- 
netic field B reveals the relative magnitude of the asso- 
ciated energy gaps, 

gap > gap „ gap , n 

e j/=o > e v=±2 e i/=±i,±3- y'- L ) 

In particular, Feldman et al, , via conductance mea- 
surements in suspended bilayer graphene, observed full 
degeneracy lifting for B > 3T, and noted that the v — 0, 
v = — 2 and v = — 1 states become apparent at 0.1 T, 0.7 



T and at 2.7 T, respectively. Similarly, in bilayer devices 
on SiC>2/Si substrates Zhao et al? 2 observed full degen- 
eracy lifting for B > 20 T and noted that the resistance 
minima for v = — 2 and v = — 3 are barely affected by 
a parallel magnetic field B\\ . Both experiments indicate 
that the u = — 1 and v = — 3 QH states emerge at similar 
magnetic fields. These observed features are consistent 
with our picture of the pseudo-zero-mode octet for case 
(i) Hz > u ~ 0. 

The energy gaps inferred or extracted from exper- 
iments 21-25 are gene rally smaller that those expected 
from V c s=s (70/eb) \J B[T] meV. In particular, Martin 
et al., 24 via local compressibility measurements of sus- 
pended bilayer graphene devices, probed energy gaps of 
size el% w 1.7 B[T] mcV and ef^ ±2 « 1.2 B[T] meV 
(and also e^2±i ~ 0.1 B[T] with less data points); also a 
recent transport experiment 25 reports a somewhat larger 
gap ef^ Q ~ 5.5 B[T] meV. Such linear B scaling of gaps 
may appear to contradict the naive \[B behavior of V c . A 
possible resolution is that charge is efficiently screened 27 
due to quantum fluctuations of the valence band; e& is 
more enhanced for smaller B. This screening weakens V c 
and could entail quasi- linear B behavior of V c . 

Let us try to handle some numbers. Supposing that 
e STo ~ 1 7 me V at 10T (using the numbers quoted above) 
and putting it into Eq. (6.7) yields V c ~ 14 meV or 
£6 ~ 15, i.e., quite a suppression of V c . On the other 
hand, according to Eqs. (6.6) and (6.7), e^ — e^_ 2 u =° 
Hz + \(z)W > Hz- Supposing e^!? - e^_ 2 ~ 5 mcV 
at B =10T then yields ~ 3. This suggests that the 
difference ef^ — e^_ 2 contains a significant Coulombic 
contribution other than the Zeeman energy (though the 
effect of charge screening is not clear). This is a mere 
order-of-magnitude estimate and one has to remember 
that experimental numbers themselves have considerable 
uncertainties; the value of tb easily changes with the in- 
put one uses. At least, the observed values of energy 
gaps indicate that they are due to the Coulomb interac- 
tion, which, though screened, still sets a far larger scale 
than the spin and valley breakings. 

The effect of an electric field E± on broken-symmetry 
states has also been studied 23 ' 25 for suspended bilayer 
graphene. Wcitz et al. 23 measured the two-terminal con- 
ductance as a function of density and electric field at var- 
ious magnetic fields. They find, in particular, that the 
v = state shows quantized conductance except near two 
large values of E±, and interpret this as suggesting the 
crossover of the spin-polarized v — state at low E± to 
layer-polarized v — states at large E±. The data also 
suggest a similar crossover at large E± for the v = ±1 
states. 

The crossover phenomena are seen in our picture as 
well. In our picture the v = state is different in com- 
position for u ~ and u ^> u CI . The crossover of the 
spin-polarized v — state to the valley-polarized state, 
when E± is applied for a fixed density, has to proceed by 
changing their composition. As illustrated in Fig. 4 (a), 
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FIG. 4: The pseudo-zero-mode sector (a) at v = and (b) at 
^ = — 1. Green arrows indicate how levels with valley (+, — ) 
move with increasing interlayer bias u. 



the filled (0+j,, 1+j.) levels go up and the empty (0_+, 1_+) 
levels come down with increasing u. A crossover thus 
takes place when the gap is closed, e? + ~ 0, i.e., around 40 

u ~ 2c 1 (0); conductance quantization is naturally lost 
around the zero gap. An observed value E±_ ~ 50 mV/nm 
for B — 2.65T, read from Ref. 23, e.g., amounts to u ~ 17 
meV, which indeed is on the order of magnitude of the 
Coulombic gap. A similar crossover around u ~ 2e^(0) is 
also expected for the v = — 1 state; see Fig. 4 (b). 



mode sector alone, has to be properly taken into account 
in studying the dynamics in the pseudo-zero-mode sector. 

Another rinding of the paper is negative capacitance 
energies, which suppress possible valley rotations in the 
pseudo-zero-mode sector. They come from the nontrivial 
structure of the interlayer Coulomb interaction in bilayer 
graphene, such that the n = 1 modes are distributed 
on both layers while the n = modes are localized on 
either layer; the valley and layer degrees of freedom are 
generally not the same for bilayer graphene. 

Experimentally it will be a challenge to study the 
structure of the pseudo-zero-mode levels directly via mea- 
surements of cyclotron resonance. Cyclotron resonance 
obeys the selection rule 41 A|n| = ±1 and is diagonal in 
spin and valley. There is no cyclotron resonance within 
the pseudo-zero-mode sector at v = and v = ±2 be- 
cause of a mismatch in spin or valley. Cyclotron reso- 
nances in the 1 +-> channel will emerge around v = ±1 
and v = ±3 with resonance energies [ranging from Ae v 
to (1 — 5cf)py considerably higher than previously ex- 
pected. 13 Such resonances are in practice rather difficult 
to detect because the pseudo-zero-modes carry current 
much smaller (by factor zM — z(u/2)/lu c <C 1) than 
other \n\ > 2 modes, as read from the Hamiltonian 



VIII. SUMMARY AND DISCUSSION 



e£oj c 
'1/5 



McAA^RZ.o+ASRZJ (8.1) 



In a magnetic field bilayer graphene has an octet of 
zero-energy Landau levels due to an extra twofold de- 
generacy in Landau orbitals n — and n = 1; this de- 
generacy is not accidental and has a topological origin. 
These levels evolve, in the presence of spin and valley 
breakings and Coulomb interactions, into pseudo-zero- 
mode Landau levels. 

In this paper we have studied the structure of this 
pseudo-zero-mode octet and, in particular, pointed out 
that its orbital degeneracy is lifted by Coulombic quan- 
tum fluctuations of the Dirac sea (the valence band). 
This quantum effect derives from the Dirac-sea contri- 
bution to the electron selfenergy, and the n = and 
n = 1 modes are split due to their different spatial dis- 
tributions. It is analogous to the Lamb shift, 33 which is 
the energy difference between the 2P 1 / 2 and 2S 1 / 2 states 
of the hydrogen atom due to vacuum fluctuations. 

For bilayer graphene vacuum fluctuations make the 
ri = modes higher in energy than the n = 1 modes 
(of the same spin and valley) when they are empty, but 
the order is reversed when they are filled. We have seen 
that such vacuum effects are intimately correlated with 
the Coulomb interaction within the zero-mode sector so 
that the spectra of this special sector correctly realize 
the particle-hole symmetry of the basic Hamiltonian. In 
a sense, a multi-component version of the Lamb shift 
emerges in bilayer graphene owing to spin and valley de- 
grees of freedom. It will be clear now that this quantum 
effect, though simply overlooked in earlier approaches 
that rely on Coulomb interactions projected to the zero- 



for a time- varying potential A — A x (t) t 
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iAy(t), where 
One can equally well ex- 



plore the structure of the pseudo-zero-mode sector by 
looking into the n = — 2 — > 1 and 1 — > 2 channels of cy- 
clotron resonance at different filling factors v £ [—4, 4]. 
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Appendix A: rotations in valley and (0, 1) space 

In this appendix we outline the derivation of the ex- 
change interaction V£ z in Eq. (4.8). Let us start with 
v x 

the kernel Q is rewritten as 



J2p v pll^- m E q- (4- 1 )- Via rotations (4.4), 



[M«)<" + M«)<1C 



l:bal 



~yab 



(Al) 



»:6t 



where Nf° a cx (G\<f>%° $™ ;a |G) denote the filling-factor 



matrices of the level n € (0,1), T^ a g.Q are the charge 



operators for the transformed fields, i.e., R a p.' with 
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-+ *f, and X% = < ' /<;':;;;' + e^Tl'X- The 
coefficient functions are defined as 
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4(\g 
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2sgCg 



+ 44 n 



Hi 2 

K°\ 2 - 



° 2 4 (c^ 2 -I c ^ + 2S 



HfliT-lfl? 



00|2\ 



for v = 0. The suppressed terms (• • •) oc eventually 
vanish when is chosen to make t(9) stationary; we omit 
such terms from now on. 

The exchange interaction AV^ Z inherits the structure 
of Vl z in Eq. (4.1). It consists of two parts: AV^\ = 

2 S p v p 7p + ^lfl-»7i) f° r tnc ^-K"' + JT'.K' combination 

and AV^ - £ p « p 7 2 (n - Q\ g ^ h ) for the tftf + K'K' 
combination, where f2| 9 _>h stands for £1 in Eq. (Al) with 
gmn re pl aC ed by ftp™. For calculations one may use 

(A2) -^l' (-^) anc ^ note *^at ^Isr^/i = ^p 1 — ffp 1 = — ^! to ~ 
gether with 



where q = £\p\, sg = sin(#/2) and eg = cos(#/2) for short. 
The special combination 5, on substituting g™ n , turns 
out to vanish. Noting the formula X) P w p7p [1, Q 2 , 9 4 ] = 
[1, 1, 3]y/n/2 V c and carrying out integrals over p with 
functions /„(<?) and f n {6) leads to V£ z in Eq. (4.8). 



Appendix B: 0(v p ) correction AV pz 

In this appendix we construct the 0(v~) ~ 0(14 d/t) 
correction AV pz to the Coulomb interaction V pz acting 
within the n = (0, 1) sector. The 0(d/t) direct interac- 
tion reads 



AV^ = -\Wv. 



nn;aa 
act 



R 



■Mb 



+ \W (/C 8v™) K m 5Rp% (Bl) 



]T « P 7p[l, 9 2 , q 4 ] = (W/2) [1, 2, 8] + 0(d 2 /0- (B3) 



The result for AV| Z to 0(W) is 



t(9)w[(F^SN° a + F^5N} ia ) 5K n a % 
+ I{F°(0)<W° Q - F^SNjJi SX afj ] , (B4) 



with ne (0,1); sn^ = n n ^ ++ 
XZj--X-f,etc. Here 

F°(6>) = 1 -D(l -cos6>) 2 
= F»{9)=c\-{D 



K 



F^(9) = 1-25(1 



cost 



1 - -2 

- sin t 
jS) sin 2 i 

- \z sin 2 i 



(B5) 



with n, m e (0, 1); = v™> KK - v™> K ' K ' and W = 

2p v~ =0 = {d/t) V c . This, via rotations (4.4), is rewritten 
as 

AV^ Z = -\vWlt^ h 

+ {l-t(§)}W{I(6) + ±K(6)} + ---,(B2) 



with D = ±c((l-c() 
Similarly one finds 



AV£ Z _ 



\~zW{t{9)h i , 



1 -t'(9)I 2 }, 



h =Ii- t{9){5N$ a - 5N} a }{5Kl% - STZ^. ), 
h =12- t(9)(SN° 0a - 8N} ja ) 6X a0 , (B6) 



with 1(9) and K{9) defined in Eq. (5.10); t{9) = 
i(sin^) 2 . The first term is proportional to the filling 
of the n — (0, 1) sector and is here adjusted to vanish 



with X\ and Xi defined in Eq. (5.10). Collecting AV^ Z , 
AV^ + and AV^_ and combining terms involving t(9) into 
capacitance energies yields Eq. (5.9). 
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